Solidity of viscous liquids. V. Long-wavelength dominance of the dynamics 
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This paper is the fifth in a series exploring the physical consequences of the solidity of highly 
viscous liquids. Paper IV proposed a model where the density field is described by a time-dependent 
Ginzburg-Landau equation of the nonconserved type with rates in k space of the form To + Dk 2 . 
If a is the average intermolecular distance, the model assumes that D ^> To a 2 . This inequality 
expresses a long-wavelength dominance of the dynamics, which implies that the Hamiltonian (free 
energy) to a good approximation may be taken to be ultralocal, i.e., with the property that equal- 
time field fluctuations are uncorrelated in space. Paper IV also briefly discussed how to generalize 
the model by including the molecular orientational fields, the stress tensor fields, and the potential 
energy density field. In the present paper it is argued that this is the simplest model consistent 
with the following three experimental facts: 1) Viscous liquids approaching the glass transition do 
not develop long-range order; 2) The glass has lower compressibility than the liquid; 3) The alpha 
process involves several decades of relaxation times shorter than the mean relaxation time. The 
paper proceeds to list six further experimental facts of viscous liquid dynamics and shows that these 
follow naturally from the model. 
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I. INTRODUCTION 



The idea that viscous liquids approaching the calori- 
metric glass transition are solid-like goes back in time 
at least to Kauzmann's and Goldstein's famous papers 
from 1948 and 1969 [HQ. According to these authors, 
when a molecule changes position in a highly viscous liq- 
uid, this happens in the form of a so-called flow event, 
a sudden rearrangement of a group of molecules. In this 
picture, which was later confirmed by computer simula- 
tions [H, 0, H, d, 0] , all motion is vibrational in the time 
between flow events. This indicates that a viscous liquid 
is much like a disordered solid p|. 

The property that highly viscous liquids are solid-like 
and more to be viewed as "solids that flow" than like or- 
dinary less- viscous liquids, is here termed solidity. This 
paper is the fifth in a series (I-IV) H, [Mill El devoted 
to extracting the physical consequences of the solidity of 
highly viscous liquids. A discussion of solidity, its motiva- 
tion and consequences, may be found in the introduction 
to paper IV to which the reader is referred for further 
physical background. 

A crucial ingredient of solidity is time-scale separation 
in the equilibrium viscous liquid. This is the fact that, 
whereas some processes occur on the vibrational time 
scale, i.e., over picoseconds, the genuine relaxation pro- 
cesses are much slower. Depending on temperature the 
latter occur on time scales of milliseconds, second, days,... 
The below discussion focuses exclusively on modelling the 
relaxation processes. 

One paper utilizing arguments from solid-state elastic- 
ity to viscous liquid dynamics preceded this series. This 
was a joint publication with Olsen and Christensen from 
1996 where the "shoving model" for the temperature de- 
pendence of the viscosity (or relaxation time) was pro- 



posed 13[. According to this and related elastic models 
the activation energy is proportional to the in- 
stantaneous shear modulus Goo (the shear modulus mea- 
sured on a very short time scale). This elastic constant 
increases significantly upon cooling, enough to explain 
the observed non-Arrhenius behavior for several molecu- 
lar liquids [H, 0, EB| (more data, however, are needed 
to illuminate whether this generally explains the non- 
Arrhenius viscosity). 

The first and second papers of this series focussed on 
the individual flow events. In paper I a "solidity" length 
I was introduced characterizing the length scale below 
which there is time between two flow events to establish 
clastic equilibrium. In terms of the average intermolecu- 
lar distance a, the alpha relaxation time r, and the high- 
frequency sound velocity Coo, the solidity length is given 
by l A = a 3 rc 00 . Close to the calorimetric glass transition 
the solidity length approaches 1 /xm, so glass-forming liq- 
uids are solid-like on quite large length scales. The model 
discussed in paper IV and here focuses on dynamics be- 
low the solidity length. 

Papers III and IV [ll|, Q3 dealt with the alpha relax- 
ation process and how to explain its seemingly generic 
high-frequency loss oc uj^ 1 / 2 16]. Paper III approached 
the problem inductively by noting that the BEL model 
from 1967 [ItJ fits data for the frequency-dependent shear 
modulus well. Starting from the BEL model it was ar- 
gued that the w -1 / 2 high-frequency behavior arises from 
a long-time-tail mechanism operating over a range of 
times shorter than the alpha relaxation time. This was 
justified by a solidity-based argument with the irrele- 
vance of momentum conservation as an important ingre- 
dient, thus allowing for the centre of mass to move follow- 
ing a flow event. As detailed in Refs. [18|, [l9( momentum 
conservation is irrelevant for highly viscous liquids, just 
as it is irrelevant in theories for defect motion in crystals. 
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Paper IV took a deductive approach and proposed a 
field-theoretic model giving a concrete realization of the 
idea of paper III. In this model, which is the subject also 
of the present paper, the generic uj^ 1 / 2 high-frequency 
decay of the alpha loss derives from a third-order term in 
the Hamiltonian (free energy) [2(|. The model, besides 
regarding momentum conservation as irrelevant just as 
in other stochastic models (e.g., for polymer dynamics 
[2lj). was based on the observation that density conser- 
vation is also apparently disobeyed. This follows from 
solidity: A flow event brings the liquid from one poten- 
tial energy minimum to another. Any potential energy 
minimum corresponds to a state of elastic equilibrium, 
thus a state of zero divergence of the stress tensor - a 
"solid" state. From the standard equations of elastic- 
ity [12] it follows that upon a local density change, the 
leading term in the displacement field in the surround- 
ings is radial and varies with distance r from the flow 
event as 1/r 2 . This is a pure shear displacement [22| . 
implying that in a coarse-grained description density can 
change at one point without changing elsewhere; den- 
sity has the appearance of not being conserved. This is 
conveniently condensed into the following equation for 
the coarse-grained density dynamics (where b^ is a di- 
mensionless measure of the magnitude of the flow event 
taking place at r M at time t^): 

p(r,t) = 536 M 5(r-r M )<y(t-t M ). (1) 

Equation |T]) does not constitute a theory, of course, since 
it does not describe how flow events correlate; it just 
describes the equilibrium coarse-grained density fluctu- 
ations. Nevertheless, Eq. (p} serves to emphasize that 
density has the appearance of not being conserved, a re- 
sult which is less trivial than the already-mentioned mo- 
mentum nonconservation deriving directly from the ex- 
tremely large kinematic viscosity [23l characterizing liq- 
uids approaching the glass transition [19( . The model dis- 
cussed in paper IV and below incorporates the main idea 
behind Eq. {!]), density nonconservation, into a frame- 
work that is explicitly consistent with statistical mechan- 
ics. The model description of the density dynamics is 
more involved than that of Eq. ([1]), which is too sim- 
ple because it corresponds to isotropic flow events. On 
very large length scales, though, the model dynamics are 
described by Eq. {J). 

In the present paper, supplementing the arguments 
of papers III and IV, we adopt a combined induc- 
tive/deductive approach. First, we list three experimen- 
tal facts characterizing highly viscous liquids from which 
the model is arrived at as the simplest model consistent 
with these facts; this constitutes the inductive part of the 
paper (Sees. 2 and 3). In Sec. 4 we proceed to discuss 
six further experimental facts and their relation to the 
model. It is argued that most of these are consequences 
of the model whereas some, though not mathematical 
consequences, appear quite natural when viewed in light 



of the model. This constitutes the deductive part of the 
paper. Section 5 gives a brief discussion. 



II. THREE FACTS OF VISCOUS LIQUID 
DYNAMICS 

This section lists three facts characterizing glass- 
forming liquids approaching the calorimetric glass tran- 
sition. 

• Fact 1: There is no long-range density- coupled or- 
der 

A popular approach to understanding why the viscosity 
increases ten orders of magnitude for a temperature de- 
crease of typically just 10-15% is to assume that some sort 
of long-ranged order gradually develops upon supercool- 
ing. According to several prominent models the dramatic 
relaxation-time increase is a consequence of the liquid ap- 
proaching a critical point where the relaxation time be- 
comes infinite (see, e.g., the excellent recent reviews by 
Tanaka (24j, Tarius and coworkers [25|, and Lubchenko 
and Wolynes [26j| . and their references). Following the 
theory of critical phenomena it is assumed that there is a 
diverging correlation length at the critical point. There 
is no consensus on how to define the proposed diverging 
correlation length, though, i.e., which quantity develops 
long-ranged correlations. Numerous X-ray and neutron 
scattering experiment however show that whatever this 
hypothetical quantity may be, it does not couple to the 
density field: No increase of long-range density fluctua- 
tions is observed upon cooling [27j |. 

• Fact 2: The glassy phase has lower compressibility 
than the liquid phase. 

The glass transition is a falling-out-of-equilibrium tak- 
ing place when the liquid relaxation time becomes 
longer than the characteristic laboratory time scale. 
Compliance-type linear-response quantities like specific 
heat, compressibility, and thermal expansion coefficient 
all decrease when going from the liquid to the glass [27| • 
This is easy to understand, because if each compliance- 
type linear-response quantity has contributions from 
both the fast (vibrational) and the much slower (configu- 
rational) degrees of freedom, these linear-response quan- 
tities must decrease at the glass transition since below 
T g the configurational degrees freeze and cease to con- 
tribute. 

• Fact 3: The alpha process is characterized by 
a distribution of relaxation times covering several 
decades of times shorter than the alpha relaxation 
time. 

Dielectric relaxation experiments are often fitted by the 
frequency-dependent response function corresponding to 
the stretched exponential dipole time-autocorrelation 
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function, exp[— (t/r)' 3 ] 28]. This reflects the fact that 
viscous liquids with few exceptions [29j do not have sim- 
ple exponential time-autocorrelation functions. There is 
nothing magic about the stretched exponential; it gives 
a good single-parameter fit to data because it reproduces 
the observed loss peaks that are always asymmetric to the 
high-frequency side [28]. The point to be made here is 
that if one rewrites any good fit to the observed autocor- 
relation function as a sum of exponentials, the distribu- 
tion function must include several decades of relaxation 
times shorter than the main (alpha) relaxation time. 



III. THE SIMPLEST MODEL FOR 
EQUILIBRIUM VISCOUS LIQUID DYNAMICS 
CONSISTENT WITH FACTS 1-3 

In this section we consider the question: Which vari- 
able^) must be included in a useful theory? What are the 
simplest dynamics for this/these variable(s) consistent 
with facts 1-3? Most models for viscous liquid dynamics 
attempt to explain both the alpha relaxation characteris- 
tics at a given temperature and the non-Arrhenius tem- 
perature dependence of the alpha relaxation time. We 
here leave aside the non-Arrhenius problem altogether, 
although solidity appears to play an important role for 
this property as well [Hj]. The idea is to "cut the Gor- 
dian Knot" by separating the complex problem of viscous 
liquid dynamics into two independent, hopefully easier 
problems. 

Nowadays most of physics - from particle physics to 
critical phenomena, electromagnetism, and condensed- 
matter physics in quite diverse contexts - is formulated 
in the language of field theory. It is natural to expect 
that viscous-liquid dynamics should also be described by 
a field theory [30(. The question which variables are rel- 
evant thus becomes: "Which fields must be included in 
the description?" 

The obvious fields are those of standard hydrodynam- 
ics: the momentum, energy, and particle density fields. 
As mentioned, the conservation laws for momentum and 
energy are both irrelevant for viscous liquids [H [lj]. 
Thus the standard hydrodynamic description, based on 
continuity equations for these quantities, looses its physi- 
cal significance. Despite the fact that solidity also implies 
(apparent) density nonconservation, the situation is dif- 
ferent for the particle density field - after all, molecules 
are not continuously exchanged with the surroundings in 
the way that momentum and energy are. We thus base 
the model sought for on the density field p(r,i). The 
question which other fields to include in order to have 
a complete description of the macroscopic dynamics, is 
dealt with at the end of this section. 

The next question is: What are the simplest pos- 
sible dynamics? Density dynamics of viscous liquids 
have two parts, the vibrations (phonons) and the "re- 
laxing" part of the dynamics. At low temperatures the 
alpha relaxation time is much larger than picoseconds, 



and the two dynamics are very well separated. Thus it 
makes good sense to ignore the vibrational part of the 
dynamics. Assuming that the relevant field theory is 
based on a Hamiltonian (i.e., a free energy functional), 
the question next is how to model the dynamics consis- 
tent with the Hamiltonian. The answer to this is well- 
known, use Langevin dynamics that is the generic way 
to arrive at dynamics from statics [2l|, [3l|: If the rel- 
evant variables are denoted by Q\, Q n , the Hamilto- 
nian is H (Qi, Q n )i an d = 1/ksT, the Langevin 
equation is Qi — —Tidi((3H)/dQi + £i(t) where is 
a Gaussian white noise term with zero mean obeying 
= 2Ti5ijS(t — t'). These equations repro- 
duce the canonical probability oc exp(— j3H\ thus ensur- 
ing consistency with statistical mechanics [U [3l| . 

The system consists of N molecules with coordinates 
Tj in volume V, and the density field is defined by 
p(r) — J2j3( r — r j)- As always when there is transla- 
tional invariance it is convenient to go to k-space; the 
range of allowed k-vectors is limited to the discrete set 
consistent with periodic boundary conditions. We define 
the k'th density component as 



(2) 



With this normalization pk fluctuations become indepen- 
dent of volume for V — > oo (N/V = Const.) and the re- 
laxational part of the static structure factor S(k) is given 
by S(k) = (pk/O-k) for V — > oo. The Langevin equation 

mlmlis 



= -r, 



(3) 



The complex Gaussian white noise term obeys ££(^) = 
£_ k (t) and (&(*)&(*')) = 2T k 6 Kk ,5(t - t'). Equation 
Q is a standard time-dependent Ginzburg-Landau equa- 
tion. Because p k = p_k this equation in conjunction with 
the equation for p_k is equivalent to two real Langevin 
equations, one for the real part of pk and one for its 
imaginary part. 

Following standard field theory procedure we split H 
into a sum of a quadratic "free-field" term Hq and an 
"interaction" term H' containing all higher-order terms: 



H 



Hn + H'. 



(4) 



Regarding H' as a perturbation, let us focus on the dy- 
namics embodied in the free-field time-autocorrelation 
function denoted by (pk(0)p-k(t))o- If there were no 
higher order terms, because S(k) = (pkP-k), the free- 
field Hamiltonian would be given by 



/3H 



S(k) 



(5) 



4 



Substituting this into Eq. ([3]) we find that the free-field 
time-autocorrelation function (pk(0)p-k(i))o 1S an expo- 
nential with decay rate 7& given by 



When H' is just a perturbation, the distribution of relax- 
ation rates is roughly given by the distribution of 7/c's. 

Since molecules cannot disappear, one would a pri- 
ori assume Tk cx k 2 ("conserved" case - model B of 
Ref. (Hi), reflecting the expectation that density at suf- 
ficiently long wavelengths obeys the diffusion equation. 
This, however, is inconsistent with fact 2 for the following 
reason. At the glass transition the liquid high-frequency 
compressibility becomes the glass compressibility. Fact 
2 states that this quantity is lower than the liquid (dc) 
compressibility, implying that in the equilibrium liquid 
there are relaxational volume fluctuations on a macro- 
scopic length scale taking place on a finite time scale (the 
alpha time scale). (This is also known from measure- 
ments of the frequency-dependent bulk modulus show- 
ing that the low-frequency bulk modulus is smaller than 
the high-frequency bulk modulus [33[, as well as from 
the so-called Mountain peak (relaxation mode) of light 
scattering [34], |35|.) This implies that the relaxational 
part of the static structure factor, S(k), has a non-zero 
limit for k -> 0. Therefore, if T fc oc fc 2 , by Eq. © the 
rate of the relaxational macroscopic density fluctuations 
would go to zero for fc — -> 0, i.e., be arbitrarily slow for 
large enough samples. This violates fact 2. Consequently, 
density fluctuations cannot be described by Tk oc k 2 as 
k — > 0, and density must have the appearance of a non- 
conserved field. 

Given that oc k does not work, the simplest alter- 
native is that Tfc is independent of k: = Tq ("non- 
conserved" case - model A of Ref. [Hj]). This, however, 
does not work for the following reason. For ka ~ 1 the 
static (relaxational) structure factor S(k) is of order one, 
whereas for k — ► S(k) converges to the ratio between 
liquid (relaxational) compressibility and that of an ideal 
gas at the same density. This ratio is typically of order 
10" 2 [H. Thus if T k = T Q , Eq. © would imply a range 
of relaxation times covering at most 2 decades. This is 
inconsistent with fact 3 that requires several decades of 
relaxation times. (It would also be physically counterin- 
tuitive to have long-wavelength density fluctuations de- 
caying much faster than short- wavelength fluctuations.) 

Since oc k 2 contradicts fact 2 and = Tq contra- 
dicts fact 3, the next possibility is a combination of the 
two: 

r fe = r + Dk 2 . (7) 

For this to be consistent with fact 3, however, D must be 
quite large: Unless the Dk 2 term makes T^ vary several 
decades for the range of allowed k-vectors, the model 



won't work for the same reason that = Tq doesn't 
work. Since the maximum k obeys ka ~ 1, this means 
that the following must be assumed: 

D > r a 2 . (8) 

The inequality ([5]) expresses a long-wavelength domi- 
nance of the dynamics because it implies that the Dk 2 
term dominates the rate expression Eq. |(7J) for a range of 
small k vectors corresponding to wavelengths much larger 
than a. Whereas in papers III and IV this inequality was 
justified by microscopic arguments, it here comes about 
from a search for the simplest possible model consistent 
with facts 1-3. 

In conjunction with fact 1, long- wavelength dominance 
of the dynamics implies that static (equal-time) correla- 
tions of density fluctuations at differing points can have 
little influence on the dynamics. Consequently, in the 
name of simplicity such correlations will be assumed to 
be absent altogether. A field theory with no equal-time 
correlations between fields at different points in space is 
termed "ultralocal." In k-space the assumption of ul- 
tralocality means that all coefficients are k independent 
in the expansion in orders of p^: Hq = ^ k (l/2A)pkP-k 

and H> = Ek,k'(V3VA0p k Pk'P-k-k' + - M- 

A similar long-wavelength dominance of the deeply 
supercooled dynamics is crucial also in Schweizer's and 
Saltzman's theory of activated hopping in polymer melts 
[HI]. Specifically, atomistic structure in S(k) is here 
coarse-grained over, resulting in a statistical segment 
level (nm-scale) description where confining forces are 
quantified via the amplitude of thermal density fluctua- 
tions, which are proportional to the isothermal compress- 
ibility, i.e., given by S(0). In other words, it is assumed 
that the most important dynamic length scale is signifi- 
cantly larger than a monomer. 

The assumption that equal-time density fluctuations 
are uncorrelated in space does not mean that the model 
is trivial - there are non-zero correlations between den- 
sity fluctuations at different positions at different times. 
The situation is analogous to that of a trivial spin model 
(i.e., with no spin-spin interactions) with Kawasaki dy- 
namics, where a spin flip at one point can take place only 
if a neighboring spin simultaneously flips in the opposite 
direction. If by chance a given up spin is surrounded by 
a large island of up spins, the given spin will be frozen for 
some time; thus the dynamics of this spin is influenced 
by those of its surrounding spins even though there are 
no equal-time spin-spin correlations. 

As a minimum, a theory for viscous liquid dy- 
namics should make it possible to calculate all 
macroscopic frequency-dependent linear-response quan- 
tities. Response functions are determined from the 
fluctuation-dissipation theorem by first calculating time- 
autocorrelation functions of variables like the total dipolc 
moment (dielectric constant), pressure (bulk modulus), 
shear stress (shear modulus), or energy density (specific 
heat). Note that the frequency-dependent bulk modulus 
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may be determined from density fluctuations alone, im- 
plying that the pressure is not an independent variable 
in the present context, and the stress tensor has only 5 
relevant components. 

These considerations lead (paper IV) to the following 
general recipe for modelling viscous liquid dynamics: 

1. The relevant degrees of freedom are fields (f>^(r), 

cf>( n >(r) defined as: a) the densities of the differ- 
ent molecules, b) the densities of the molecules' var- 
ious configurational variables reflecting the molec- 
ular symmetry, c) the 5 stress tensor fields of the 
traceless stress tensor, d) the potential energy den- 
sity field; 

2. The Hamiltonian H (free energy) is ultralocal and 
consists of scalar (i.e., rotationally and translation- 
ally invariant) terms; 

3. For each field the dynamics are described by a 
time-dependent Ginzburg-Landau equation, i.e., a 
Langevin equation of the form 



Ai) 



(*) 



(9) 



where £^ (t) is a Gaussian white-noise term obey- 



ing 



(*)) 



f) 



4. For each density field the Langevin equation 



rates are given by T 
->0')„2. 



DU') k 2 where 



£)0') ^> Tq a 2 ; for all other fields the rates are 
k-independent: =T^\ 

Inclusion of the extra fields makes it possible to cal- 
culate the frequency-dependent dielectric constant and 
shear modulus, as well as all 8 fundamental frequency- 
dependent thermoviscoelastic response functions [391 ] : 
the isochoric and isobaric specific heats, the isothermal 
and adiabatic compressibilities, the isobaric expansion 
coefficient, the adiabatic contraction coefficient, and the 
isochoric and adiabatic pressure coefficients. 



IV. SIX FURTHER FACTS OF VISCOUS 
LIQUID DYNAMICS AND THEIR 
INTERPRETATION IN LIGHT OF THE MODEL 

Below are listed six experimental facts that are each 
nontrivial in the sense that there is no logically com- 
pelling connection between it and facts 1-3. Neverthe- 
less, when viewed in light of the model, these facts all 
appear obvious. 

• Fact 4: Below the alpha-loss peak the loss virtually 
follows the Debye prediction, i.e., is proportional to 
UJ. Thus there is an effective cut-off at long times 
in the relaxation-time distribution. 



Alpha loss peaks, which are typically measured by di- 
electric relaxation experiments, a re g enerally asymmetric 
towards the high-frequency side [28| . At low frequencies 
they follow the Debye prediction [ex 1/(1 + iuir)), i.e., 
the loss is virtually proportional to frequency; the same 
applies, e.g., for the frequency-dependent shear modulus 
[40|. This implies that, if the linear response is written 
formally as a sum of Debye processes, there is an effective 
cut-off at long times in the relaxation-time distribution. 

Given the spatial disorder of any viscous liquid it is 
not surprising that Debye peaks are rarely observed, but 
the observed effective cut-off at long relaxation times is 
highly nontrivial. If the non-Debye relaxation were due 
to effects of disorder-induced activation-energy broaden- 
ing, the obvious guess would be a Gaussian activation- 
energy distribution. This, however, implies loss peaks 
that are symmetric in a log- log plot, which is inconsis- 
tent with experiment. 

Because there is a minimum relaxation rate in the 
model, a long relaxation-time cut-off is automatically im- 
plied. Of course, this cut-off was built into the model via 
Eq. ([7]). The point we wish make is merely that this 
equation was not justified from fact 4, but as the sim- 
plest way to rationalize facts 1-3. Fact 4 follows. 

• Fact 5: Above the alpha loss-peak frequency the 
loss appears to be generically close to lo^ 1 / 2 . 

The alpha process is conveniently monitored by dielec- 
tric relaxation experiments [28]. Above the loss-peak 
frequency the dielectric loss decays following an approx- 
imate power law oc oj~ v '. The shape of the loss peak 
often changes with temperature, i.e., time-temperature 
superposition (TTS) is often not obeyed. It is now 
generally agreed that Johari-Goldstein beta processes 
may be found at much lower frequencies than previously 
thought [4l[. Thus, since alpha and beta relaxations 
have quite different temperature dependence, low-lying 
beta processes easily lead to TTS violations even in the 
Hertz regime and below. According to this reasoning the 
"generic" characteristics of the alpha process is observed 
only when TTS is accurately obeyed. A study of the di- 
electric loss of 10 molecular liquids published in 2001 [l6| 
indicated that when TTS is accurately obeyed, the expo- 
nent v is quite close to 1/2. A recent study comprising 
45 molecular liquids found that the minimum slope a m i n 
above the loss peak (characterizing the inflection point 
and thus giving the best approximate power-law fit) is 
generally fairly close to —1/2 (60% of the liquids studied 
obey -0.6 < a min < -0.4) g|. 

As shown in Ref. [2(| and paper IV, to lowest order 
in perturbation theory a third-order term in the Hamil- 
tonian implies that the loss varies as uj~ x I 2 for oj's con- 
siderably higher than the loss peak frequency. 

• Fact 6: The alpha relaxation process is dominated 
by small-angle jumps. 

Many models assume "cooperatively rearranging re- 
gions" that do not interact with one another. In this 
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spirit molecules either move significantly (those involved 
in the flow event) or do not move at all (those in the sur- 
roundings). Thus molecular jumps would be expected 
to be fairly large and molecular orientations likewise to 
change considerably following a flow event. The molec- 
ular jump angles cannot be probed by linear-response 
experiments, but fortunately they can by NMR experi- 
ments. The result of Bohmer and coworkers [43l |44| is 
that small-angle jumps dominate. 

This observation calls for an explanation in terms of so- 
lidity: The reason that small-angle jumps dominate must 
be that - if the flow event picture is not completely wrong 
- the overall picture is dominated by the small adjust- 
ments in the surroundings required to re-establish elastic 
equilibrium after a flow event. A simple solidity-based 
calculation (paper I) shows that the jump-angle distribu- 
tion P(4>) varies as 1/0 2 , consistent with NMR findings 
[43l ] (this distribution is not normalizable because there 
are infinitely many molecules in the surroundings - in 
reality the distribution is cut-off at very small angles be- 
cause elastic effects do not propagate beyond the solidity 
length). Clearly, largest weight is given to small jump 
angles. 

When fact 6 is contemplated in light of the model, 
it should be noted (paper IV) that the crucial To 7^ 
identity expressing apparent density nonconservation can 
come about only if a flow event is followed by small 
solidity-based adjustments of molecular positions in the 
far surroundings. Thus fact 6 follows from the model's 
nonzero IV 

• Fact 7: Viscous slowing down is not accompanied 
by significant changes of the static structure factor. 

A popular and obvious explanation of the dramatically 
increasing relaxation time for liquids approaching the 
calorimetric glass transition is that upon cooling there 
is a gradual build-up of some sort of long-range order, in 
many models signalling that there is a critical point not 
far below T g where the relaxation time becomes infinite 
HH [26[ . As mentioned, numerous experiments have 
looked for long-range order, but found none [13] ■ Not 
only is there no long-range order, but the liquid structure 
as probed by S(k) via X-ray or neutron scattering exper- 
iments changes little over the temperature range where 
the relaxation time changes by 10 or more orders of mag- 
nitude. It is matter of taste whether or not one regards 
this as surprising [451 ] - A model with no spatial correla- 
tions at any temperature trivially predicts that there are 
no nontrivial changes of the static structure factor upon 
cooling. In this sense fact 7 follows from the model. 

• Fact 8: The Debye-Stokes-Einstein relation is of- 
ten violated in viscous liquids. 

An important finding of the 1990's was that translational 
motions often decou ple from and become 1-3 decades 
faster than rotations [461 ] . Somehow, translations are en- 
hanced compared to what one expects from the Debye- 
Stokes-Einstein relation that estimates the single particle 



diffusion constant D s from the viscosity 77. Although this 
relation a priori applies only for macroscopic particles, 
there is no obvious reason that a molecule on average 
should move much longer than a during the rotational 
correlation time (the alpha relaxation time). 

The generally accepted picture of Debye-Stokes- 
Einstein violations is that these reflect dynamic hetero- 
geneity [461 ] . Although dynamic heterogeneity is not de- 
scribed by the model, this is not inconsistent with the 
model that provides only a coarse-grained description 
of the dynamics. Debye-Stokes-Einstein violations, in 
fact, fit nicely with the model if one assumes that the 
single-particle diffusion constant is roughly equal to the 
diffusion constant D of Eq. ([7]); in this case the long- 
wavelength dominance inequality ([5]) simply expresses 
Debye-Stokes-Einstein violation. 

• Fact 9: The dynamics are not sensitive to the 
chemistry. 

At any given temperature the relaxation characteristics 
- linear as well as nonlinear - are similar for all glass- 
forming liquids. Viscous liquids do differ as regards, e.g., 
whether or not there is a clearly defined beta relaxation 
in the liquid phase, or whether or not time-temperature 
superposition applies. But these differences do not ap- 
pear to correlate with chemical characteristics in any sim- 
ple fashion. For instance, as regards their dielectric and 
thermodynamics properties, liquids held together exclu- 
sively by van der Waals forces do not differ systematically 
from those involving hydrogen bonds [28| . Similarly, the 
frequency-dependent shear modulus of viscous metallic 
liquids forming bulk metallic glasses is indistinguishable 
from that of typical molecular liquids [47[. The over- 
all conclusion is that it is not possible from macroscopic 
measurements alone to determine which chemical bonds 
are involved (except, of course, from the fact that the 
temperature range where the glass transition takes place 
trivially provides information about the strength of the 
intermolecular forces). As a further confirmation of fact 
9 it should be noted that mixtures behave much like pure 
liquids and that, again, it is not generally possible from 
purely macroscopic measurements to tell whether or not 
a given liquid is a mixture. 

This insensitivity to the chemistry is a highly signif- 
icant fact. For some reason the dynamics are fairly 
indifferent to details of the molecular interactions, but 
why? In the model chemistry independence is expected 
because of the long- wavelength dominance of the dynam- 
ics - clearly chemistry plays little role for dynamics on 
length scales much larger than the size of a molecule. 
The chemistry independence is similar to that of critical 
phenomena. This has motivated many attempts to draw 
parallels to the theory of critical phenomena by assuming 
that viscous slowing down is accompanied by some kind 
of long-range order responsible for the observed quasi- 
universality. The present model has no such assumption, 
but on the contrary assumes that there are only short- 
range static correlations of the relevant fields. Only when 
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it comes to the dynamics do long length scales play im- 
portant roles. - Note that the model does not imply ab- 
solute chemistry independence, because the free energy 
as function of the fields varies with chemical composition. 

V. DISCUSSION 

It is a long-standing assumption that viscous- liquid dy- 
namics are cooperative. Mode-coupling theory [43, |49| 
is an interesting case where cooperativity enters via the 
coupling of single-particle motion to the surroundings, 
resulting in a modification of the single-particle motion 
with drastic consequences at low temperatures (infinite 
relaxation time at a finite temperature in the simplest 
and most studied version of the theory). The present 
model is also cooperative, but in a rather simple-minded 
sense, with two elements of cooperativity: 1) Density 
nonconservation (implying Tq ^ 0) is a cooperative ef- 
fect because it involves the solidity-based movement of 
molecules far from a flow event; 2) The long-wavelength 
assumption Eq. ([8|) expresses cooperativity in the sense 



that motion over long wavelengths is mainly responsible 
for the alpha process. 

The overall purpose of this paper was to show that 
the model of paper IV has implications that were not 
put into the model. The overall credibility of the model 
is strengthened by the fact that the model is consistent 
with - and in most cases predicts - several experimental 
facts that are independent of the model inputs. We would 
like finally to emphasize that, despite the macroscopic 
reasoning of this paper, the long-wavelength dominance 
inequality Eq ([8]) reflects properties of the individual 
flow events [50 ]. 
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